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Abstract 

We present a theoretical study of the motion of the antihydrogen atom {H) in the Earth's gravi- 
tational field above a material surface. We predict that H atom, falling in the Earth's gravitational 
field above a material surface, would settle in long-living quantum states. We point out a method 
of measuring the difference in energy of H in such states that allow us to apply spectroscopy 
of gravitational levels based on atom-interferometric principles. We analyze a general feasibility 
to perform experiments of this kind. We point out that such experiments provide a method of 
measuring the gravitational force (Mg) acting on H and they might be of interest in a context of 
testing the Weak Equivalence Principle for antimatter. 

PACS numbers: 
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I. INTRODUCTION 



Galileo, Newton and Einstein recognized that all bodies, regardless of their mass and com- 
position, fall towards the Earth with an equal gravitational acceleration. Is that conclusion 
valid for antimatter? This has never been tested. 

In the context of the general relativity theory, the universality of free fall is often referred 
to as the Weak Equivalence Principle (WEP). Violations of WEP could occur in ordinary 
matter-matter interactions e.g. as a result of the difference in the gravitational coupling to 
the rest mass and that to the binding energy. WEP is being tested with increasing sensitivity 
for macroscopic bodies. The best test so far confirms WEP to the accuracy of 2-10^^'^ (using 



a rotating torsion balance Ongoing projects aim at the accuracy of 1 part in 10^^ (laser 
tracking of a pair of test bodies in a freely falling rocket [2]), or even to 1 part in 10^* (in 
an Earth orbiting satellite [3j]). However, in view of difficulties in unifying the quantum 
mechanics with the theory of gravity, it is of great interest to investigate the gravitational 
properties of quantum mechanical objects, such as elementary particles or atoms. Such 
experiments have been already performed, e^g. using interferometric methods to measure 



the gravitational acceleration of neutrons ^,l5| and atoms 



with antz atoms ( see 
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However, the experiments 



11| and references therein) are even more interesting in view of 



testing WEP, because the theories striving to unify gravity and quantum mechanics (such 
as supersymmetric string theories) tend to suggest violation of the gravitational equivalence 
of particles and antiparticles 12| . Experiments testing gravitational properties of antiatoms 



are on th e ag enda of all experimental groups working with antihydrogen (see e.g. ATHENA- 
ij, ATRAP Q and AEGIS [jj). he challenging aspects in experiments 



ALPHA 



of this kind is to control the initial parameters of antiatoms, such as their temperature and 
position, with sufficient accuracy |16 |. 

In the present paper we investigate a possibility to explore gravitational properties of 
antiatoms in the ultimate quantum limit. We study antihydrogen atoms levitating in the 
lowest gravitational states above a material surface. The existence of such gravitational 



17H19[|. The existence of analogous states 



states for neutrons was proven experimentally 
for antiatoms seems, at a first glance, impossible because of annihilation of antiatoms in the 
material walls. However, we have shown that ultracold antihydrogen atoms are efficiently 
reflected from material surface 20|, l2l|] due to so-called quantum reflection from the Casimir- 
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Polder atom-surface interaction potential. We have shown that antihydrogen atoms, confined 
by the quantum reflection via Casimir forces from below, and by the gravitational force 
from above, would form metastable gravitational quantum states. They would bounce on a 
surface for a flnite life-time (of the order of 0.1 s) [20]. This simple system can be considered 
as a microscopic laboratory for testing the gravitational interaction under extremely well 
specifled (in fact, quantized!) conditions. 

The annihilation of ultraslow antiatoms in a wall occurs with a small but flnite (few 
percent) probability. It provides a clear and easy-to detect signal, which might be used 
to measure continuously the antiatom density in the gravitational states as a function of 
time. If antiatoms are settled in a superposition of gravitational states, the antiatom density 
evolves with beatings, determined by the energy difference between the gravitational levels. 
The transition frequencies between the gravitational levels are related to the strength of 
the gravitational force Mg, acting on antiatoms; here M is the gravitational mass of H, 
and g is the Earth's local gravitational fleld strength. Also we show that a measurement of 
differences between the energy levels would allow us to disentangle Mg in a way independent 
on effects of the antiatom-surface interaction. 

The plan of the paper is the following. In section [TTl we study the main properties of 
the quasi- stationary gravitational states; in section IIIII we present the time evolution of 
the H gravitational states superposition; in section |IV] we discuss a concept of a quantum 
ballistic experiment, namely the spatial-temporal evolution of the H gravitational states 
superposition, in section |V] we analyze the feasibility of measuring H atom properties in 
gravitational quantum states. In the Appendix we derive useful analytical expressions for 
the quasi-stationary gravitational states scalar product. 



II. H GRAVITATIONAL STATES 

In this section we discuss the properties of the H gravitational states above a material 
surface. 

We consider an H atom bouncing above a material surface in the Earth's gravitational 
fleld. H is conflned due to the quantum scattering from the Casimir-Polder potential below, 
and the gravitational fleld above. The Schrodinger equation for the H wave-function "^{z) 
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in such a superposition of atom-surface and gravitational potentials is: 



2'mdz'^ 



+ V{z) + Mgz - E 



^{z)=0 (1) 




Here z is the distance between the surface and the H atom, and V{z) is the atom-surface 
interaction potential with a long-range asymptotic form V{z) ~ —C^/z"^. We distinguish 
between the gravitational mass, that we refer to as M and the inertial mass, denoted by m 
hereafter. The wave-function "^{z) satisfies the full absorption boundary condition at the 
surface {z = 0) [21I], which stands for the annihilation of antiatoms in the material wall. 
The characteristic length and energy scales are 

(2) 
(3) 

Here Iq = 5.871 fim is the characteristic gravitational length scale, Icp = 0.027 fim 
is the characteristic Casimir-Polder interaction length scale, Eq = 2.211 10~^^ a.u. is the 
characteristic gravitational energy scale, and Ecp = 1.007 10^^ a.u. is the Casimir-Polder 
energy scale. As one can see, the gravitational length scale is much larger than the Casimir- 
Polder length scale Iq ^ Icp, while the gravitational energy scale is much smaller than the 
Casimir-Polder energy scale Eq <^ Ecp- It is useful to introduce the gravitational time scale 

To = h/Eo ~ 0.001s (6) 
For large atom-surface separation distances z ^ Icp the solution of eq.([T]) has a form: 

^W~Ai(f --) + ir(E)Bi(f --) (7) 

where Ai(x) and Bi(x) are the Airy functions [22]. The requirement of square integrability 
of the wave-function — 00) — )• results in the following equation for the energy levels 
of the gravitational states in the presence of the Casimir-Polder interaction: 

K{En) = (8) 
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n 








1 


2.338 


1.407 


13.726 


2 


4.088 


2.461 


24.001 


3 


5.521 


3.324 


32.414 


4 


6.787 


4.086 


39.846 


5 


7.944 


4.782 


46.639 


6 


9.023 


5.431 


52.974 


7 


10.040 


6.044 


58.945 



TABLE I: The eigenvalues, gravitational energies and classical turning points of a quantum bouncer 
with a mass of (anti) hydrogen in the Earth's gravitational field. 

The hierarchy of the Casimir- Polder and gravitational scales Icp suggests that the 
quantum reflection from the Casimir- Polder potential can be accounted for by modifying the 
boundary condition for the quantum bouncer (a particle bouncing in the gravitational field 
above a surface, the interaction of the latter with a particle is modeled by infinite reflecting 
wall). The quantum bouncer wave-function satisfies the following equation system: 



2mdz'' 



+ Mgz - En ^n{z) = 

^ ^n{z = 0) = 

The quantum bouncer energy levels are known to be equal to 



(9) 



hi 



El 



Ai(-A° 



0. 



(10) 

(11) 



Table [T] summarizes the eigenvalues and classical turning points z° = E^/ (Mg) for the 
first seven gravitational states of a quantum bouncer (with the mass of antihydrogen). 

For the distances Icp ^ z Iq one could neglect the gravitational potential in eq.([T]). 
In this approximation, the solution of eq.([T]) has the following asymptotic form: 



^(z) ~ sm{kz + 5{E)). 



(12) 



Here k is the wave vector k = \/2mE, and S{E) is the phase-shift of H reflected from the 
Casimir- Polder potential in absence of the gravitational field Matching asymptotics 
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in eq. (fT2|) and eq.([7]) we get a relation between the phase-shift S{E) and the i^— function 
introduced in Eq.(I7]): 

^ tan((5(E)) Ai-(-i^/eo) - kloAi{-E/eo) 
^ ' tan((5(E))Bi'(-E/£o) - A;/oBi(-E/£o)' ^ ' 

In deriving the above expression we took into account that the relation between K{E) and 
S{E) should not depend on the matching point and thus can be formally attributed to 
Zm = 0. An equation for the distorted gravitational levels could be obtained by substitution 
of eq.(Il3]) into eq.®: 

tan{S{En)) ^ Ai(-En/eo) 

kh M'{-Ejs,y ^ ' 

This equation is equivalent to the following boundary condition: 

$(0) _ tan(^(K)) 

Thus the following equation system describes an H atom, confined by the Earth's gravita- 
tional field and the quantum reflection from the Casimir-Polder potential: 



2md. 



-:^ + Mgz-En M^) = , , 

-I (16) 

■^■(0) ^ tan(5(£;„)) ^ ' 
$'(0) k 



For the lowest gravitational states the condition klcp <^ 1 is valid. Thus the scattering 
length approximation for the phase-shift ~ —kacp is well justified. The complex-value 
quantity [21i] : 

acp = -(0.10 + il.05)/cp, (17) 
acp = -0.0027 - i0.0287/im (18) 

is the scattering length on the Casimir-Polder potential provided full absorbtion in the 
material wall. 

Thus the equation for the lowest eigenvalues iHM has a form: 

ocp ^ Ai{-En/eo) , . 

lo Ai'i-EjeoY ^ ^ 

The above equation is equivalent to the following boundary condition for the wave- 
function $(2) of a particle in the gravitational potential eq.(IT]): 

<!>{z ^ 0) ^ z - acp (20) 



Because of the imaginary part of the scattering length acp-, the gravitational states of H 
above a material surface are quasi-stationary decaying states. For low quantum numbers n, it 
is easy to relate the lowest quasi-stationary energy levels to the unperturbed gravitational 
energy levels of a quantum bouncer. Indeed the variable substitution z = 1 + acp 
transforms Eq.( 19pn]) to the equation system for the quantum bouncer: 



+ Mgz-{Er,-Mgacp) 



= (21) 



2md'z^ 

<l'„(?^0)-^0 (22) 

The eigenvalues En and eigenfunctions $„ are : 

En = El + MgacP, (23) 
$„(^) = i-Ai((^-acp)//o-A°), (24) 

where Ni is the normalization coefficient (see Eqs. ( ITQllHOj) in the Appendix). In the following, 
we will use the dimensionless eigenvalues A„ = En/sQ. 



\n = Xl + acp/lo (25) 

An important message from the above expression is that the complex shift Mgacp (due 
to the account of quantum reflection on the Casimir-Polder potential) is the same for all low- 
lying quasi- stationary gravitational levels. It means that the transition frequencies between 
the gravitational states are not affected by the Casimir-Polder interaction, provided the latter 
can be described by the complex scattering length acp- The scattering approximation is 



valid in the limit knttcp — 0, where kn = \/2mE^ (let us note that for the first gravitational 
state |fciacp| = 0.0071). However, accounting for the higher order fc-dependent terms in 
Eq.f ll4l) would result in the state dependent shift of the gravitational states due to the 
Casimir-Polder interaction. We use a known low energy expansion of the s-wave phase-shift 



S{E) in a homogeneous l/z"^ potential 23|, in which we keep the two leading ^-dependent 
terms: 

dcp{k) cot {6{k)) ^ -1 + ^^{Icpk) + Ulcpkf In (^-^) + ... (26) 

6 acp o \ 4 / 

We introduce a k-dependent modified "scattering length" acp{k) = —6{k)/k and get the 
following expression for acp{k): 

acpik) ~ acp + '^Icpilcpk) + ^acpilcpkf In (27) 
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The leading k-dependent term in the above expression ^lcpi}cpk) is real and independent 
on properties of the inner part of the Casimir-Polder interaction. It is determined by the 
asymptotic form of the potential, thus it depends on the Casimir-Polder length scale Icp 
only. Then the modified equation for the gravitational state energies is: 

En = El + Mg~acp{En) (28) 

Taking into account the smallness of the k-dependent terms (for the lowest gravitational 
states) in expression (1271) . we get: 

E^c^EU Mg~acp{K) =e[\l + acp/k + ^(^cpfc°) + ^(^cpA:°)^ In ^) . (29) 

Here = ^J2mE^. 

The account for k-dependent terms in Eg. (1271) modifies the transition frequencies be- 
tween the gravitational states in a way, dependent on the Casimir-Polder interaction. How- 
ever, such modification is very weak. Indeed, taking into account, that Icpk^ ~ Icp/Iq 
for the lowest gravitational states , the leading k-dependent term corrections to the grav- 
itational energy are of the second order in a small parameter Icp/Iq- The transition fre- 
quency between the first and second gravitational states equals uji2 = cj?2 + ^12? where 
cjO^ = (^0 - E^)/{2nh) = 254.54 Hz, and A12 = Mg{dcp{k^2) - acp{k^i)) = 0.0017 Hz. 

An account of the first two terms in Eq. (!27|) provides equal decay width for the lowest 
gravitational states. This width is determined by the probability of antihydrogen penetrating 
to the surface and annihilating. 

r. = e^. (30) 

Here we use a standard notation b = 4Imacp: 

b = 0.115 (xm. (31) 

The widths of the gravitational states fl30|) are proportional to the ratio s/Iq. Using Eqs. ([2]) 
and (jH) we could find that this ratio is equal to the gravitational force s/Iq = Mg so that 

r„ = ^-Mg. (32) 
The corresponding life-time (calculated for an ideal conducting surface) is 
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r=-^~0.1s. (33) 
Mgb ^ ' 

We note factorization of the gravitational effect (appearing in the above formula via a 

factor Mg) and the quantum reflection effect, manifestating through the constant h. Such 

a factorization is a consequence of the smallness of the ratio of the characteristic scales 

h/{2lo) ~ 0.01. 

Comparing the H lifetime in the lowest gravitational states with the classical period 
T = ~ 0.0033 s of H with the ground state energy bouncing in the gravitational 

field, we see that H bounces in average about 30 times before annihilating. This shows that 
the lowest gravitational states are well resolved quasi-stationary states. 

It is interesting to estimate the maximum gravitational quantum number A^, below which 
the gravitational states are still resolved, i.e: 



Tn r(iv)^ 



> 1. (34) 

dn 

Here tn is the lifetime of the A^— th gravitational state. Tat is a classical period, corresponding 
to the A^— th state via T/v = 2h7T /{dE{n)/dn). For such an estimation we transform eq. ffT^ 
using the asymptotic form of the Airy function for a large negative argument and get: 

\ 2/3 

2V.~ 4.-'^(K))J • (35) 

The accuracy of the above equation increases with increasing n; it gives the energy value 
within a few percent even for n = 1. In the energy domain of interest |5(-E')| ^ 7r[n — j] 



2M, so: 

An - All - (36) 



Here we used the semiclassical approximation for 24j : 

A°-(|vr[^-j) (37) 



One could verify that in the case of small n the above equation reduces to eq. (12511 . Substi- 
tution of eq. (l36|l into eq.( 1M|l results to: 

— ^ f-ik) 

T„-4Im5(E0)- ^ ^ 
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FIG. 1: The number of H bounces during the hfetime of n— th gravitational state 

The ratio T{n)/T{n) expresses the number of H classical bounces during the lifetime of the 
n-th state. This dependence is shown in Fig{TJ Using numerical values S{E), calculated in 



2l| . we find that inequality ( l3ll) holds for 

n<N = 30000. (39) 

The corresponding energy = 6 10~^^ a.u., and the characteristic size of such states is as 
large as Hn = 1.6 cm. This means that the concept of the quasi-bound gravitational states 
is justified not only for the lowest states, but it also might be applied for highly excited 
states. 

The quasi- stationary character of the antiatom gravitational states above a material 
surface manifests itself in a nonzero current through the bottom surface {z = 0). Indeed, 
the expression for the current is: 

.(^,0 -*•(.. (40) 

taken at 2; = for a given gravitational state turns out to be equal to: 

Here A„ is given by Eq. (!25l) . Ni is the normalization factor. We take into account the 
smallness of the ratio acp/lo and Eq.( ITT]) . and get : 
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Ai(-A„,)^-^Ai'(-A°), (42) 

which is exact up to the second order in the ratio acp/lo- Now taking into account an 
exphcit form of the normahzation coefficients (Eq.l lHOll in Appendix) Ni = Ai'((— A°), we get 
finally: 

This result is in full agreement with Eq.f pOj) as far as: 

jj^ mz,t)\'dz=j{0,t) = --exp{--t). (44) 

III. BOUNCING ANTIHYDROGEN 

In this section, we are interested in the evolution of an initially prepared arbitrary su- 
perposition of several lowest gravitational states of H. In the following, we will limit our 
treatment to the scattering length approximation for describing the Casimir-Polder interac- 
tion, and will neglect all, except the first, term in the expression (127|) . so that acpik) ~ acp- 
The corresponding H wave-function is: 

'^(^'^) = E^^K^Ao - A,)exp(-a,-). (45) 

Here tq is the characteristic H bouncer time scale, Ci are expansion coefficients and Ni = 
Ai'(— Ai) are the normalization factors of the gravitational states (see the Appendix). 

We are interested in the evolution of the number of antihydrogen atoms as a function of 
time: 

POO " POO (J*(J. 

F{t) = / \^z,t)\Hz=y] / -^M*{z/k-\,)ki{z/k-\i)eM-H\-^])t)dz 
Jo ij^^ Jo 

(46) 

First, let us note that the above expression for the total number of particles is no longer 
constant because of the decay of the quasi- stationary gravitation states. Second, the quasi- 
stationary gravitational states corresponding to different energies are non-orthogonal: 

1 f°° 
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In the Appendix we will derive the following expression for the cross-terms aij, exact up to 
the second order of the small ratio acp/lo- 

V(2/o) 

= ^0-AO + .6/(2/o) ^''^ 

As one can see, such cross-terms vanish if there is no decay, i.e. if 6 = Almacp 0. 

Now we can calculate an expression for the number of antihydrogen atoms as a function 
of time (Uni): 



Fit) ^ exp(-£.) {± ICf + 2 Re I E e.p(-.(A» - A»)l) 

(49) 

From Eqs. fHOj) and ( 130|) we get the following expression for the disappearance (annihila- 
tion) rate — keeping the terms up to the second order in the ratio acp/W 



\ i i>j j / 

For a superposition of the two gravitational states with the equal coefficients Ci,2 (say 
Ci = C2 = 1), the above expression gets a simple form: 

= -I exp(-^t) (1 + cos{u^2t)) . (51) 

Here UJ12 = (A2 — A'j')/ro. The same result could be obtained by calculating the flux j(0,t) 
Eq.( HO|) for a superposition of states ( H5|) . 

One can see, that the disappearance rate decays as a function of time according to the 
exponential law with the width F (the same for the lowest states), also it oscillates with 
the transition frequency between the first and second gravitational states (equal to 254.54 
Hz). We plot in Figj2] the time evolution of H disappearance rate in a superposition of 
two lowest states. Curiously, the oscillation of disappearance rate is the direct consequence 
of decaying character of gravitational states. Indeed, such an oscillation is observable due 
to the nonvanishing contribution of the interference term in the expression for the total 
probability to find antihydrogen atoms, given by Eq.(l49l). As one can see from Eq.( l48l) this 
contribution is proportional to the imaginary part of the scattering length, it would vanish 
in case there were no decay of gravitational states due to annihilation in the material wall, 
described by parameter b/lo- 
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FIG. 2: Evolution of the annihilation rate of H atom in a superposition of the first and the second 
gravitational states. 
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FIG. 3: Evolution of the annihilation rate of H atom in a superposition of first, second and third 
gravitational states. 

So far the oscillation frequency of the disappearance rate N{t) corresponds to the energy 
difference between the unperturbed gravitational levels. Expression ( l50l) does not include 
the shift of gravitational state energies Re acp/lo as it is equal for all the gravitational states, 
thus it is canceled out in the energy difference. The account for higher order k-dependent 
terms in fl27|) would result in a small (second order of the ratio (acp/^o)) correction to the 
transition frequency. A measurement of the oscillation frequency uJi2 given by Eq. (ISTll 
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would allow us to extract the following combination of the gravitational and the inertial 
masses from eq.(jl]): 

Under the additional assumption of the equality of the known inertial mass of the hydrogen 
atom niH and that of antihydrogen, imposed by CPT, we get: 



The evolution of three gravitational state superposition provides information not only 
about the characteristic energy scale Eq but also about the level spacing as a function of 
quantum number n, characterized by the value d'^E{n) / dn^ . Such a study rnight be in- 
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26| and 



teresting for testing additional (to Newtonian gravitation) interactions (see 
references there in) between H and a material surface with the characteristic spatial scale 
of the order of micrometers. Such interactions would manifest as nonlinear additions to the 
gravitational potential, which would modify the spectrum character. In the case of three 
state superposition, the disappearance rate (l50l) has the form: 



dFi23(t) 2r r /3 , , , , \\\ 

— ■ exp(-— t) - + COS(Wi2t) + COS(W23^) + COS((Wi2 + U)23jtj ■ (54) 



dt 3h h ' \2 

Here Uij = (A° — \^)/tq. One could verify that the period of coherence of cos(u;i2t) and 
cos(ci;23t) terms is: 



27r , , 

Tr = ~ 0.02s. (55) 

Ui2 — UJ23 



A semiclassical expression for Tr is: 



WEjd^y ^^^^ 



One can see that the period is a quantum limit analog of a half revival period Trev = 
4TT/\d'^E/dn'^\ {Trev characterizes the time period after which the evolution of the wave- 
packet returns to the semiclassical behavior, see ^] for details and reference therein). In 
Figj3]we plot the annihilation events as a function of time (|5^ for a superposition of three 
lowest gravitational states. The period Tr is clearly seen as a period of modulation of a 
rapidly oscillating function. The ratio 
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is sensitive to any nonlinear addition to the gravitational potential. Indeed while linear 
corrections to gravitational potential can only change eo, nonlinear additions change the 
derivative of levels density \(PE/dn^\. 



IV. QUANTUM BALLISTIC EXPERIMENT 

Two independent experiments are needed in order to determine the gravitational mass 
M and the inertial mass m of antihydrogen. In the previous section, we showed that a 
combination of gravitational and inertial masses M and m, given by Eq. (l52|) . can be ex- 
tracted from the frequency measurement eq.( l5Tl 1541) . An independent information could be 
obtained from measurement of the spatial density distribution of in a superposition of the 
gravitation states, for instance, in the flow-throw experiment (a kind of a beam scattering 
experiment), in which ff-atoms with a wide horizontal velocity distribution move along the 
mirror surface. The time of flight along the mirror should be measured simultaneously with 
the spatial density distribution in a position-sensitive detector, placed at the mirror exit. 
Such a detector would be able to measure the density distribution along the vertical 
given time instant. The horizontal component of H motion could be treated classically. Due 
to a broad distribution of horizontal velocities in the beam, atoms would be detected within 
a wide range of time intervals between their entrance to the mirror and their detection at 
the exit. In such an approach, we could study the time evolution of H probability density 
at a given position z: 

|$(i2)(^,t)r = exp(-^t) (|<l>^[2)(^)l' + 2Re$;^*)(^)exp(-za;i2t)) (58) 



h 

ki{z/k - Ai) 



Ai'(-Ai; 



Ai(z//o - A2) 



Ai'(-A2) 



(59) 



^int _ Ai(^Ao-Ai)Ai(2://o-A2) , . 

^^''^^'^ - Ai'(-Ai)Ai'(-A2) ^ ^ 



The transition uji2 = 254.5411 z could be extracted from the probability density time 
evolution at a given z. The length scale /q could be extracted from the position of the zero 
z\ of the wave-function. Here superscript stands for the state number, and a subscript 
corresponds to the number of zero for a given state. Thus such a position is determined by 
the condition Ai{zi /Iq — A2) = 0; is equal to the following expression: 

) = (A2 - Ai)/o = 10.27/im. (61) 
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FIG. 4: Color. The probability density of -ff in a superposition of the first and second gravitational 
states, as a function of the height z above a mirror (vertical-axis) and the time t (horizontal-axis). 
Dark shade: low probability density. Light shade: high probability density. The dashed line 
indicates the position of the node in the wave-function of the second state 

The probability density in Eq.f l58|) of a two states superposition ok, z = z\ behaves like the 
probability density of the ground state alone: 



The probability density at a height z\ does not exhibit any time-dependent oscillations. 
We show the probability density as a function of the height z above a mirror (y-axis) and 
the time t (x-axis) in a superposition of the first and second gravitational states in FigHl 

The position of the node in the wave-function of the second state is shown in FigJH as 
a horizontal line, separating lower and upper rows of periodic maxima and minima in the 
probability density plot. The position z\ does not depend on initial populations of the 
gravitational states, which makes it beneficial for extracting the spatial scale /q. 

From knowing the length Zq and the energy eo scales, one could get the following expres- 
sions for the inertial m and gravitational M masses of H: 




(62) 



m 



260^0 



2 ' 



(63) 



M 



£o_ 



(64) 
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The equality m = M postulated by EP relates Eq and Iq as follows: 

^0 = ^yj^' (65) 
or, using the gravitational time scale Eq.(l6]): 

.^/f. (66) 

One can easily recognize in the above expression a classical time of fall from the height Iq 
in the Earth's gravitational field. 

Thus a measurement of the temporal-spatial probability density dependence of in a 
superposition of the two lowest gravitation states would provide a full information on the 
gravitational properties of antimatter. The superposition of three (and more) gravitational 
states could be useful to search for additional (to gravity) interactions with a spatial scale 
of the order of Iq- For such a purpose, it is useful to study the probability density at zeros 
of each Airy function in the state superposition. In particular, the positions corresponding 
to zeros of second and third gravitational states are the following: 

= (A2 - Ai)/o = 10.27/im, (67) 
zf) = (A3 - A2)/o = 8.41/im, (68) 
= _^ ^(3) ^ i8_Q8^rn, (69) 

The three state (ijk) superposition probability density at position of zero is equal to the 
two state superposition (jk) probability: 

l%.)(4'=\t)p = |$(,,)(zf,t)p. (70) 

This means that the three state probability density exhibits harmonic time-dependent oscil- 
lation with a frequency Uij at height of zero z^. Let us mention that the time dependence 
of |$jjfc(z,t)p in any position z, except for the mentioned zeros is not harmonic; it is given 
by a superposition of three cosine functions with different frequencies, analogous to Eq llMll . 
This property allows us to extract the zeros positions using the probability density. One can 
see that the knowledge of zeros (Eq.( |68ll69l) ) is analogous to the knowledge of the transition 
frequencies. A measurement of one zero position allows us to extract the spatial scale Iq, a 
measurement of positions of two or more zeros allows us to constrain hypothetical nonlinear 
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additions to the gravitational potential. We show the probability density as a function of 
height z above the mirror (y-axis) and the time t (x-axis) in a superposition of first, second 
and third gravitational state in FigJ^i 




5 10 15 

t(ms) 



FIG. 5: Color. The probability density of in a superposition of the first, second and third 
gravitational states, as a function of the height z above a mirror (vertical-axis) and the time t 
(horizontal-axis). Dark shade: low probability density. Light shade: high probability density. The 
dashed lines indicate the positions of the nodes in the wave-functions of the second and third state. 



V. FEASIBILITY OF ANTIHYDROGEN GRAVITATIONAL STATES EXPERI- 
MENT 

In this section we study the feasibility of an experiment on the gravitationally bound 
quantum states of anti-hydrogen atoms. For such an estimation, we compare it with the 
already performed gravitational states experiments using ultra-cold neutrons (UCN) 



19j. The UCN gravitational experiments can be used as a benchmark for such a comparison 
because 1) the neutron mass is nearly equal to the anti-hydrogen mass, 2) any modification 
to the H quantum state energies and wave-functions following from the precise shape of the 
Casimir-Polder potential are small compared to those of the quantum bouncer; UCN in the 
Earth's gravitational field above a perfect mirror is well described by the quantum bouncer 
model, 3) our estimation of B. lifetimes in the quantum states of 0.1 s are compatible to or 
even longer than the time of UCN passage through the mirror-absorber installation, also 4) 
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UCN velocities are comparable to velocities of ultra-cold H atoms produced in traps Q ■ 
We will discuss here mainly the statistical limitations arising from an estimate of the spectra 
from sources of H atoms that are projected in the near future. 

In the simplest configuration, the experimental method for observation of the neutron 
gravitational states consisted in measuring the UCN flux through a slit between the hor- 
izontal mirror and the flat absorber (scatterer) placed above it at a variable height as a 
function of the slit height (the integral measuring method), or analyzing the spatial UCN 
density distribution behind the horizontal bottom mirror exit (the differential measuring 
method) using position-sensitive neutron detectors. The slit height can be changed and 
precisely measured. The absorber acts selectively on the gravitational states, namely the 
states with a spatial size = A°/o smaller than the absorber height H are weakly affected, 
while the states with if„ > H are intensively absorbed 
of the experimental method, the experimental setup, t 
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28 



29| . A detailed description 



this phenomenon could be found, for instance in 17 



le results of various applications of 
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3Q|. 



Leaving aside numerous methodical difficulties in the experiments of this kind (as they 
have been already overcome in the neutron experiments) and a real challenge to get high 
phase-space densities of trapped anti-hydrogen atoms (they are aimed at anyway in the 



existing anti-hydrogen projects 10|,lll(|), let us compare relevant phase-space densities in the 



two problems, keeping in mind that it is the principle parameter, which defines population 
of quantum states in accordance with the Liouville theorem. If the phase-space densities 
of anti-hydrogen atoms would be equal to those of UCN, we would just propose to use an 



32| for anti-hydrogen experiments with minor 



existing UCN gravitational spectrometer [31, 
modifications. 

UCN is an extremely narrow initial fraction in a much broader, and hotter neutron 
velocity distribution. Maximum UCN fluxes available today for experiments in a flow-though 
mode are equal to 4 10^ UCN / cm? / s] such UCN populate uniformly the phase-space up 
to the so-called critical velocity of about Q m/s (UCN with smaller velocity are totally 
reflected from surface under any incidence angle; thus they could be stored in closed traps 
and transported using UCN guides). If one uses pulsed mode with a duty cycle, say 10^^, 
the average flux would drop to 4 UCN/cm^/s. The pulse method provides more precise 
measurements, it is used in current experiments with the GRANIT spectrometer, and it will 
be used in gravitational interference measurements analogous to those performed with the 
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centrifugal quantum states of neutrons [33|, [3^. Taking into account the phase-space volume 



available for UCN in the gravitationally bound quantum states in the GRANIT spectrometer 



32|, we estimate the total count rate of about 10^ events/day when the relative accuracy for 
the gravitational mass is 10^'^ (we note that the accuracy in the mentioned experiment is 
defined by a width of a quantum transition, and a few events might be sufficient to observe 
the corresponding resonance). 



The average flux of H atoms projected by AEGIS collaboration [15|, |35|] is a few atoms 
per second; let's take it equal 3 H/s to have it defined. The cloud length is ~ 8 mm, its 
radius is ~ 1.5 mm, thus the cloud volume is ~ 5 x 10^^ cm?. For comparison, estimate of 
the average UCN fiux, which would be emitted from a small UCN source with a volume of 
~ 5 X 10~^ cm^, the maximum UCN density available of 30 UCN /cm?, with a duty cycle 
of ~ 10^^ gives 0.15 UCN/s. It is 20 times lower than the H fiux estimated above. One 
should not forget, however, that the projected H temperature is ~ 100 mK, i.e. 100 times 
larger than the effective UCN temperature. Thus we lose a factor of a/100 = 10 because of 
larger spread of H vertical velocities. No geometrical factors are taken into account here as 
well as no constraints following from final solid angles allowed, final sizes of H detectors, 
mirrors etc. However, their account would decrease our estimation by only a few times 
provided proper experiment design (note that equal acceleration of all anti-hydrogen atoms 
would not decrease their phase-space density). Thus we could provide statistical power of H 
experiment compatible to that with UCN. As the projected temperature of anti-hydrogen 
atoms in another proposal j^Gj is significantly lower (~ 1 mK that is just equal to the effective 
UCN temperature), the phase-space density of anti-hydrogen atoms could be even higher. 
Another significant advantage of a lower temperature consists in a more compact setup 



design. Note that a gravitational spectrometer analogous to [31|, 132| selects just a very small 
fraction of UCN( H) available (those with extremely small vertical velocity components) 
thus the count rate of "useful" events is extremely low in both cases. 

Thus, we conclude that measurements of the gravitationally bound quantum states of 
anti-hydrogen atoms look realistic if they would profit from methodical developments avail- 
able in neutron experiments plus high phase-space densities of anti-hydrogen atoms aimed 
at in future. Based on extensive analysis of the mentioned neutron experiments, we could 
conclude that measurements of the gravitational mass of anti-hydrogen atoms with an ac- 
curacy of at least 10~^ is realistic, provided that projected high H phase-space density is 
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achieved. 



VI. CONCLUSIONS 

We argue on existence of long-living quasi-stationary states of H above a material sur- 
face in the gravitational field of Earth. A typical lifetime of such states above an ideally 
conducting plane surface is r ~ 0.1 s. Quasi- stationary character of such states is due to 
the quantum reflection of ultra-cold (anti) atoms from the Casimir- Polder (anti) atom-surface 
potential. The relatively long hfe-time is due to the smallness of the ratio of the charac- 
teristic spatial antiatom-surface interaction scale Icp and the spatial gravitational scale Iq. 
We show that the spectrum of H decaying gravitational levels is quasi-discrete even for the 
highly excited states as long as their quantum number n ^ 30000. We argue that low lying 
gravitational states provide an interesting tool for studying the gravitational properties of 
antimatter, in particular for testing the equivalence between the gravitational and inertial 
masses of H. We show that, by counting the number of H annihilation events on the sur- 
face, both the transition frequencies between the gravitational energy levels, as well as the 
spatial density distribution of gravitational states superposition can be measured. An im- 
portant observation in this context is that a modification of the above mentioned properties 
of gravitational states due to the interaction with a surface disappears in the first order of 
the small ratio Icp/Iq- Finally, we show that actual measurements of quantum properties 
of H atoms, levitating above a material mirror in gravitational states are feasible, provided 
that the projected high phase-space density H is achieved. 
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VIII. APPENDIX 

Here we derive an expression for the scalar product of two complex energy gravitational 
states eigenfunctions: 

1 f°° 

aij = — — / Ai*(z//o - Xj) ki{z/k - \)dz (71) 

with 

Nf= ki^{z/h-\j)dz. (72) 
Jo 

We start with the equations for the eigenfunction A\{z/lo — Aj) and the complex eigenvalues 
-^j = -^i* + o^cp/W- 

- Ai"{z/k - \) + zAi{z/k - \) = Xi Ai{z/lo - X^). (73) 
The equations for the complex conjugated eigenfunction and the eigenvalue are: 

- Ai*"(;2//o - A,) + zAi*iz/k - Xj) = X* M*{z/lo - X^). (74) 

We multiply both sides of equation Eq.(|73ll by Ai*(z//o — A-,) and integrate them over z. 
Then we multiply both sides of Eq.( 17^ by Ai^z/Iq — Aj) and integrate them over z. After 
substraction of the results of these operations we get: 

Ai*'(-A,)Ai(-A,)-Ai*(-A,)Ai'(-A,) = (A*-A,) / ki*{z/k-X,) M{z/lo-X,)dz. (75) 

Jo 

To get the above result, we integrated by parts the integrals with second derivatives and took 
into account that Airy functions vanish at infinity. Now we take into account the equality 
Ai(— A°) = and smallness of the ratio acp/lo, to get the following expressions, exact up to 
the second order in acp/lo'. 

Ai(-A.) = -^Ai'(-A°), (76) 
Ai*(-A,) = -^Ai*'(-A;). (77) 
Up to the second order in acp/h we get: 

Ai*'(-A,) Ai(-A.) - Ai*(-A,) Ai'(-A.) = zA Ai*'(-A;) Ai'(-A°). (78) 
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The above result should be combined with the known expression for the normalization 
coefficient: 

7V2 = Ai'2(-A,) + A,Ai2(-A,), (79) 
which, up to the second order in acp/h turns to be: 



Ar. = Ai'(-A°)(l + A°|f). (80) 



Keeping first order terms, we finally get: 
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